Mon. Not. R. Astron. Soc. 000, ITIH51 (20091 Printed 8 September 2009 (MN WpX. style file v2.2) 



The UV-optical colours of brightest cluster galaxies in optically and 
X-ray selected clusters 
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ABSTRACT 

Many brightest cluster galaxies (BCGs) at the centers of X-ray selected clusters exhibit clear 
evidence for recent star formation. However, studies of BCGs in optically-selected clusters 
show that star formation is not enhanced when compared to control samples of non-BCGs 
of similar stellar mass. Here we analyze a sample of 113 BCGs in low redshift (z < 0.1), 
optically-selected clusters, a matched control sample of non-BCGs, and a smaller sample of 
BCGs in X-ray selected clusters. We convolve the SDSS images of the BCGs to match the 
resolution of the GALEX data and we measure UV-optical colours in their inner and outer 
regions. We find that optically-selected BCGs exhibit smaller scatter in optical colours and 
redder inner NUV - r colours than the control galaxies, indicating that they are a homogenous 
population with very little ongoing star formation. The BCGs in the X-ray selected cluster 
sample span a similar range in optical colours, but have bluer NUV - r colours. Among X- 
ray selected BCGs, those located in clusters with central cooling times of less than 1 Gyr are 
significantly bluer than those located in clusters where the central gas cooling times are long. 
Our main conclusion is that the location of a galaxy at the centre of its halo is not sufficient to 
determine whether or not it is currently forming stars. One must also have information about 
the thermodynamic state of the gas in the core of the halo. 

Key words: BCG; cooling flow; star formation 



1 INTRODUCTION 

As their name already suggests, brightest cluster galaxies (BCG) 
are special in that they are among the most luminous and mas- 
sive galaxies in the Universe and are located near the minimum 
of the gravitational potential of clusters and groups. Observa- 
tions show that although most BCGs are elliptical galaxies they 
are different from ordinary galaxies of that type: they have sur- 
fac e brightness profiles b etter described by two-component mod- 
els dGonzalez et al.l20"03T ). and their Fundamental-Plane projections 
have different slopes compared to normal ellipticals. This indi- 
cates the formation histori e s of BCGs are likely to be different 
dvon der Linden et alj|2007l ; iBernardi et alj[2007l : iDesroches et al.l 
l2007l ; lLiuetal.ll2008h . 

N-body simulations of BCG formation in a ACDM cosmol- 
ogy show that these galaxies form their stars early ( z~5), but as- 
semble their final masses very late ( z~0.5) <Lucia et alj J2007t) 
and references therein). The dominant stellar popultions of BCGs 
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are inde ed observed to b e old. Studies of t he ultraviolet lumi- 
nosities llHicks et alJl2005t|Pipino et all)2008l), nebular line emis- 
sion l lAllenlll995l : ICrawford et al.ll 19991: Edward et alj|2007t) . opti- 
cal absorption line spectra i Cardiel et alj 19951) . and infrared fluxes 
jEgami et alJl2006l : IO'Dea et alj|2008l) of BCGs show that a small 
amount of residual star formation is still taking pla ce in some of 
these objects. The existence of UV -brigh t cores teildfell et alj 
l2007l ; lRaffertv et al.l2008l ; |Pipino et alj2008h indicates that the star 
formation is centrally concentrated. In some cases, the star forma- 
tion and the radio jets coincid e, possibly implying that the star- 
bursts are triggered by the jets l lDonahue et alj|2007al :lo"' Pea et alj 

Eooi. 



Star formation in BCGs is often associated with the so- 
called "cooling flow" phenomenon. Star forming BCGs are often 
found nearer to the X-ray centers of clusters than quiescent B CGs 
jBildfell et ai1l2007l : ICrawford et al.ll999l : lEdward et alj2007l) and 
star formation is correlated with the cooling timescale of the gas. 
Studies using X-ray data found that BCGs with star formation 
are clearly marked by an entropy threshold of ~ 30 KeV cm" 2 
(equivalent to a cooling time scale of ~1 Gyr and indicative of a 
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strong cooling flow dRaffertv et"ai]|2008l : Icavagnolo et alj|2008al : 
IVoit et alj|2008h ). In general, the mass of young stars in the BCG 
is only a few percent of the gas mass that is predicted to be con- 
densing from the intracluster medium and accreting onto the BCG. 
This may reflect the ability of clust ers to distribute energy released 
by an AGN to the surrounding gas dVoit et al.ll2008l) . 

Around 70% of clusters selected from flux-limited X-ray sur- 
veys are classified as co oling flow clusters, and this high fraction 
holds to at least z - 0.4 fedgelll997l : |Peres et aPl 19981 : Bauer et all 
120051) . If this fraction is universal for all clusters, and enhanced 
star formation is directly related to the cooling flow phenomenon, 
then BCGs should form stars more actively than ordinary mas- 
sive elliptical galaxies. However, studies of BCGs in optically se- 
lected groups and clus ters have found that BCGs are as old as 
their field counterpa rts jCrawfordetalJl999l :l von der Linden et al. 
l2007l:lEdward et alj|2007l) . 

The aim of our study is to clarify whether BCGs are indeed 
forming stars more actively than other massive elliptical galaxies. 
We will make use of a combination of UV and optical imaging 
data. UV/optical colors are sensitive to small amounts of star for- 
mation on timescales of <100 Myr. The colors can be measured 
directly from images and will therefore not be affected by aperture 
effects. In contrast, the spectral properties of BCGs selected from 
the SDSS survey are measured through a 3 arcsecond diamteter 
fibre aperture, and hence only probe star formation in the central 
regions of the galaxies. The imaging data also allows us to study 
the colour gradients of the BCGs. We describe our sample selection 
and our measurements of 2-zone colors in Section 2. We present the 
comparison of 2-zone colors and color gradients for different BCG 
samples in Section 3. Finally, our results are discussed in Section 4 
and our conclusions are given in Section 5. Throughout this paper, 
we assu me a cosmology with H -70 km s~' Mpc~', Q.,„=0.3, and 
n A =0.7 ( dTegmark et al]|2004l) ). 



2.2 "C4 Cluster Catalog" and lvon der Linden et alj J2007t) 
catalog 

iMiller et al.H2005h developed the "C4" algorithm to identify galaxy 
clusters from the SDSS spectroscopic sample. They based their 
method on the fact that the cores of galaxy clusters and groups are 
dominated by red galaxies and used the SDSS colors(t( - g,g - r,r - 
i,i - z) for their selection. The "C4" catalog is 90% complete for 
galaxy clusters with masses greater than 2xl0 14 /h Mq ( or veloc- 
ity dispersion a greater than 500 km/s). The completeness declines 
for lower mass clusters, reaching a value of ~55% at a mass of 
10 14 /h Mq (logo- v =2.6) and -30 % at2xlO 13 /h M ( logo-,=2.4) 
( Here, we use Equation 2 from iBiviano et al. I j2006l) is used to 
convert cluster mass to velocity dispersion ). 

Von der Linden et al. (2007) (vdL07) based their study on the 
"C4 Cluster Catalog" of the Third Data Release (DR3) of SDSS 
(748 clusters in total) and carefully identified the BCGs of every 
cluster. They searched for BCG candidates within the virial radius 
and identified the BCG as the brightest galaxy closest to the centre 
of the cluster potential . Their resulting sample consisted of 625 
BCGs in groups and clusters at z <0.1. They also recalculated the 
velocity dispersions and redshifts of each cluster. 

Our study is based on the vdL07 BCG sample and we refer 
to BCGs selected from this sample as "optical BCGs". 402 of the 
vdL07 BCGs have SDSS spectra. We match the vdL07 sample with 
the fourth General Release of GALEX (GR4) and find 113 galaxies 
with MIS coverage in the NUV band (Sample 1, SI). We prefer to 
use the MIS images rather than the shallower AIS images because 
the UV emission in most of the BCGs is weak. 

As we will discuss, the majority of clusters in the vdL07 sam- 
ple are low in mass and have smaller velocity dispersions and X-ray 
luminosities than typical clusters in X-ray samples (Figure [8). We 
thus created an additional sample of 60 massive optically-selected 
clusters with log cr,. greater than 2.8 and GALEX AIS images. We 
will use this sample when we compare optically-selected clusters 
with X-ray selected clusters. 



2 DATA AND SAMPLE 
2.1 SDSS and GALEX 

The Sloan Digital Sky Survey (SDSS) l lYork et alj[2000l) has ob- 
served a quarter of the extragalactic sky and provides images in 
five photometric bands (u',g',r',i',z')> as well as spectra of galaxies 
selected from the imaging. The spectra are taken with 3" diame- 
ter fibers and cover a wavelength range from 3800 to 9100 A. The 
SDSS images have a pixel scale of 0."396 and a mean PSF (point 
spread function) of 1. "4 (FWHM). 

The Galaxy Evolution Explorer (GALEX, iMartin et al.1 

j2005h ) is an orbiting space telescope providing imaging in two 
bands: the far-ultraviolet (FUV) centered at 1528A and the near- 
ultraviolet (NUV) centered at 2271 A. The images have l."5 pixels 
and the average resolution is 4. "3 FWHM for FUV and 5. "3 for 
NUV. GALEX is performing a number of surveys distinguished 
mainly by coverage and depth. In this work we use the data from 
the All-sky Imaging (AIS) with a typical depth of 20.5 in AB mag- 
nitude and the Medium Imaging survey (MIS) with a typical depth 
of 23.5 AB magnitude. We note that a few NUV image tiles are not 
covered by FUV observations, but all FUV images have accompa- 
nying NUV images. Because the shape of the PSF near the border 
of GALEX images is distorted, we only use those BCGs within 
1200 pixels of the corresponding image center. 



2.3 Control sample 

We constructed a control sample of field galaxies matched to the 
SI BCGs, following the steps outlined in vdL07. We first sort the 
BCGs in mass and then we search the SDSS DR4 and GALEX 
GR3 cross-matched catalogues for non-BCGs that differ in redshift, 
stellar mass and g — r color by less than 0.02, 0. 1 dex and 0. 15 mag 
respectively. Each control galaxy enters the sample only once. 

Figure [T] presents a comparison of the UV and optical prop- 
erties of SI BCGs and the control sample. Note that a significant 
fraction of the highest mass BCGs lack comparable control galax- 
ies; this is because almost all the very most massive galaxies in the 
local Universe are themselves BCGs. We note that because the C4 
cluster sample is less complete at low masses, some of the control 
galaxies for these systems may themselves be BCGs. 

2.4 X-ray BCGs from lRaffertv et alj 120081) 

X-ray data provides information about the thermodynamic state of 
the gas in clusters: by deprojecting the X-ray spectra extracted in el- 
liptical annuli from X-ray images, one can derive gas temperatures 
and densities and thereby estimate gas cooling time-scales and en- 
tropy. Traditionally, those clusters with "central" (determined by 
the limiting resolution of images) cooling time-scales shorter than 
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Figure 1. A comparison of the SI BCGs with the control sample. We plot histograms of stellar mass, redshift, g — r colour, galaxy velocity dispersion, R; so 23 
radius, absolute r-band (AB) magnitude, and NUV and FUV apparent magnitudes. BCGs are in red and control galaxies are in black. In the top left-hand 
corner of each panel we list the logarithm of the Kolmogorov-Smirnov probability that the two distributions are drawn from an identical parent population (a 
99 percent probability that the distributions are different will have a value of -2). 



the age of the u niverse (~ 10 Gyr) have been classified as cool- 
ing flow clusters (Fabian 1994 ). Recently, howe ver, several studies 



dCavagnolo et al.l 



2008dlbl: IRaffertv et alj|2008l C08, R08 respec- 



tively afterwards) have attempted to measure entropy and cooling 
times at smaller cluster-centric radii using higher resolution data 
from the Chandra satellite. These studies established a more strin- 
gent cooling time threshold of ~1 Gyr (or 0.8 Gyr at a radius of 
12 kpc from IRaffertv et al.l d2008h ) below which BCGs were ob- 
served to have significantly enhanced star formation. In our study, 
we adopt the latter definition, and refer to clusters (and their asso- 
ciated BCGs) with central cooling times of less than 1 Gyr as cool- 
core clusters/BCGs, and those with central cooling times greater 
than 1 Gyr as non-cool-core clusters/BCGs, respectively. 

To assess whether X-ray and optically selected BCGs differ , 
we make use of the sample of 46 BCGs in IRaffertv et alj d2008h 



with X-ray data from Chandra. We search for optical images in 
SDSS DR7 and find 21 matches (1 1 of these have spectra ): 12 out 
of the 2 1 are in cool-core clusters and 9 are in non-cool-core clus- 
ters. In order to compare colors\stellar populations to the optically 
selected BCGs, we need to further restrict the sample to lie at z<0.1 
and have GALEX coverage, resulting in a total of only 9 clusters 
(3 non-cool-core and 6 cool-core) (see TableQJ. 



3 PHOTOMETRY 

To ensure that our photometric measurements are consistent across 
the different bands, we transform all our images to the same geom- 
etry and effective resolution so that our photometric measurements 
trace the same part of each galaxy at different wavelengths. 
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Figure 2. Some examples of typical BCGs in our dataset. The first row shows a typical BCG in our sample with low S/N in the UV, the second row a UV 
bright BCG with red colour and flat NUV - r gradient, the third row a UV-bright BCG with a steep colour gradient, and the fourth row a BCG with a blue 
NUV - r core. 



3.1 Registering the images 

We first run SExtractor 1 on each imag e to get the celestial and pixel 
coordinates for sources in the images jBertin &Arnoutsll 19961) . We 



then select the most compact sources in the NUV images and match 
them to the celestial coordinates of sources in the SDSS and FUV 
images. We run the GEOMAP/GEOTRAN 2 tasks in IRAF, with the 



http ://sextractor. sourcef orge .net/ 
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pixel coordinate pairs as input to register every image to the frame 
geometry of the corresponding NUV image. 



3.2 Convolving the images to a common PSF 

We first need to obtain the PSF for each image. For SDSS images, 
we use the "read-PS F" software 3 and read out the PSF mod- 
els for the five bands directly from the corresponding psField file. 
We transform these SDSS PSF models to a GALEX pixel scale of 
l."5. In order to get the PSF model for the NUV and FUV im- 
ages, we first select stars within a radius of 1200 pixels around the 
center of each image. After masking nearby sources, we fit a two- 
dimensional gaussian function to the stars. We cut out small stamps 
around each star and shift them so that every star is centered ex- 
actly on the center of the stamp. We add these stamps together and 
adopt the combined stack as the best PSF model for the frame. We 
check the PSF models by subtracting the stack from each of the in- 
put stars. Then we run the IMMATCH/PSFMATCH task in IRAF 
to obtain kernel functions that we use to convolve the higher res- 
olution SDSS images to the resolution of the corresponding NUV 
images. Figure [2] shows examples of the registered and convolved 
images of 4 different SI BCGs with different UV/optical colours. 
The resulting imag es are adaptively binned in two dimensions us- 
ing the algorithm of lCappellari et al. | J2003|) so that all pixels have a 
signa-to-noise ratio (S /N) above a certain fixed value. The resulting 
images are then able to show the extended, low surface brightness 
regions of the sources. The color profiles are obtained by averaging 
the colors in a set of elliptical rings with increasing radius around 
the BCG center. The rings are chosen to that all all points along 
the profile have errors below a certain fixed value. The first row in 
Figure [2] shows a typical BCG in our sample with low S/N in the 
UV, the second row shows a UV-bright BCG with red colors and a 
flat NUV - r profile, the third row shows a UV-bright BCG with a 
steep color profile, and the fourth line a BCG with a blue NUV - r 
core. 

In Figure[3] we plot the positional differences between the cen- 
troid of the BCG measured from the convolved SDSS and GALEX 
images. As can be seen, the centroids generally agree to better than 
1", but in a few cases the offsets can be as large as 3 arcseconds. We 
have checked these images by eye and we find the largest offsets 
between the centroids can be caused by (1) SExtractor centroid- 
ing on regions of star formation, which are more prominent in UV 
light, (2) failure of SExtractor to deblend the BCG from nearby UV 
bright sources, or (3) low S /N of the UV image. So it is important 
to use SDSS image decided position parameters for both GALEX 
UV and SDSS photometry to ensure consistency. 



3.3 Two-zone photometric measurements 

In order to assess whether the BCGs in our sample have colour gra- 
dients, we measure colours both inside and outside a fixed aperture 
for each object in our sample. Because the SDSS r band images 
are of good quality and high S/N ratio, we use them as the ba- 
sis for our two-zone photometry. We first mask all sources except 
the BCG. We use SExtractor to de termine the ellipticity an d posi- 
tion angle from the r band image teertin &Arnoutslll996j) . These 
determine the shape and position angle of the apertures that we 
use. Following vdL07, we define the total magnitude of the BCG 
as the light contained within the radius where the r-band surface 
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Figure 3. Differences between UV and R-band centroid positions (in 
GALEX pixel units) of the BCGs in our sample. 



brightness equals (23 + 101og(l + z)) mag arcsec~ 2 (which we de- 
note as r im 2i). All our photometry measurements are processed in- 
side this aperture, which encloses most of the light. The reason 
why we measure the BCG magnitude within r iso 23 and do not at- 
tempt to measure a total magnitude, is because the surface bright- 
ness profiles of BCGs are complex and do not follow a simple 
de Vaucouleurs profile dGonzalez et al J[2005t) . We define R 50 and 
R90 as the radii enclosing 50% and 90% of the total flux of the 
galaxy evaluated within r, !o2 3 . We measure the total UV and r-band 
flux inside R 50 (in), and between R 50 and R 90 (out). We compute 
colours within these two regions and and define the color differ- 
ence, A(NUV - r) = (NUV - r) m - (NUV - r) outs . We apply this 
definition to our UV/optical colour measurements, because it has 
higher 5 /N than more traditional measures of colour gradient. We 
caution that in Sect. 13. 51 we will use a different measure of colour 
gradient, because we are comparing directly with existing studies 
published in the literature. 



3.4 Accuracy of our 2-zone color measurements 

Recall that our two-zone photometry is carried out on optical im- 
ages that have been convolved to the resolution of the GALEX data. 
It is thus important to understand to what extent our colour differ- 
ence measurements reflect the intrinsic values for the real galaxy. 

To assess this, we simulated 10,000 galaxies in the r and NUV 
bands and included a similar level of background noise as in our 
real data (e.g. the MIS background for the simulated NUV images). 
The models assume Sersic profiles with n=4 for the r band and n=2 
for the UV band (i.e. we assume that the optical light comes from 
a bulge-like component and the UV from a star-forming disk). The 
range in r magnitude, NUV magnitude and r-band effective radius 
for the simulated galaxies are constrained to be the same as in the 
vdL07 sample. By varying the effective radius of the NUV light, 
we obtain a set of simulated galaxies with different colour gradi- 
ents. We then apply our algorithm to these images and compare our 
colour measurements to the intrinsic "input" values in Figure|4] The 
inner color is defined to be the color inside R50 and the outer color 
is defined as the colour between R 50 and R 90 - 

We begin by exploring the effect of noise of our colour mea- 
surements. The right column in Figure [4] shows the differences be- 
tween the measured colours for model images with and without 
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Figure 4. The results of tests using simulated images showing how the convolution procedure and sky noise will affect our colour measurements. The left 
column shows differences between the intrinsic (input) colours and the colours measured from the simulated images, which have been convolved to the 
resolution of the GALEX data and which have Poisson noise added. The right column shows the effect of just adding the Poisson noise. The data are divided 
into three different color bins ([5,5. 5], [5. 5,6], [6,6.5]) ( blue, green and black lines in the plot). Solid lines show the median values at a given r-band magnitude, 
while dashed lines show the 25% and 75% percentiles of the distribution at the same magnitude (see text for details). 
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Figure 5. An analysis of the influence of structural parameter differences between BCGs and non-BCGs on our 2-zone color measurements (see text for 
details). Results are shown for simulated BCGs (red) and simulated non-BCGs (black). The simulated BCGs and non-BCGs have the same distribution of 
colours. The colour profiles of both kinds of objects are assumed to be flat. The differences between the red and black histograms arise purely because of the 
structural differences between the two classes of object. 



Poisson noise. The median value of the differences is zero across 
most of the r-magnitude range spanned by our BCG sample. Sys- 
tematic effects only become important at r magnitudes fainter than 
16.1 (most of our BCGs are brighter than this). We then investi- 
gate what effect the convolution to GALEX resolution has on our 
measurements. The left column of the plot shows the measure- 
ment differences between model images with perfect resolution and 
the GALEX NUV resolution. The top panel shows that the inner 
NUV - r colour will be underestimated and the situation is worse 
for redder objects. The middle panel shows that the outer NUV - r 
colour will be over-estimated and the situation is again worse for 
redder objects. However, the outer color is better constrained in 



general than the inner colour. At an r-band magnitude of ~16.0, 
the inner NUV - r of the reddest objects will be underestimated by 
~0.3 mag, while that of the bluest ones would be underestimated by 
~ 0.1 mag. The outer NUV - r colour would be over-estimated by 
~0.15 mag on average. The bottom panel shows that the resulting 
color difference A(NUV- r) will be over-estimated by ~0.5 mag at 
r ~16.0. These results are easily understood from the fact that the 
optical light (with Sersic model n=4) is more concentrated than the 
UV light (with Sersic model n=2), and will thus be more severely 
scattered to the outer region of the galaxy. As we will show later 
on, our results are not substantially affected by these systematics. 
Von der Linden et al.(2007) found that BCGs have larger ef- 
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fective radii and lower surface brightnesses than non-BCGs. We 
simulated two sets of galaxies with structural parameters similar 
to the optical BCGs and the control sample to check the extent to 
which differences in structural parameters affect our measurements 
of 2-zone colors. We fix the r band magnitude to be 15 (~ median 
value of the optical sample), NUV - r to be 6.0 (a typical value 
at the red end of our sample) and A(NUV - r) to be 0. Only the 
mean surface brightness inside R^o, /J50, is allowed to vary (R50 will 
vary accordingly since the total r magnitude is fixed). The median 
values of yu 50 are fixed to be 18.82 and 18.65 mag for BCGs and 
non-BCGs respectively (corresponding to the median values of the 
two samples as measured by vdL07). /j 50 is allowed to vary ran- 
domly around the median value with a scatter of 0.7 mag arcsec -2 
(again according to von der Linden's distribution). We compare 
the distributions of the values of (NUV - r) m , (NUV - r) out , and 
A(NUV - r) (defined as in SectO for BCGs (red) and non-BCGs 
(black) in Figure [5] It can be seen that structural differences cause 
both the inner and outer colours of BCGs to be shift slightly blue- 
wards relative to non-BCGs. However, the color differences for the 
two samples are similar. We have also simulated cases where the 
input galaxies have an intrinsic colour gradient, and the results are 
similar. 



3.5 u-r gradients 

The SDSS optical photometry is of higher resolution than the 
GALEX photometry, so u - r colour gradients offer an alterna- 
tive way to trace star formation in the central regions of early type 
galaxies. As was found by R08 and will be seen in Sect |4.3| X-ray 
selected cool-core BCGs have steeper u-r colour gradients than 
non-cool-core BCGs. We will use u-r gradients to link our study of 
optical selected BCGs to previous studies of X-ray selected BCGs 
(e.g. R08). 

We use the SDSS r and u band images to measure the u - 
r colour profiles. These images have lower resolution and poorer 
seeing than the images used by R08, but we produce high quality 
color profiles by registering and convolving the u band images to 
the r band images. We then measure the color gradients (G(u - rj) 
following the procedures oulines in in R08. 

The procedure we employ is the following. We measure 
colours in a series of elliptical annuli, with the major axis increas- 
ing in ~ 1 pixel intervals. If the u - r profile has a positive color 
gradients near the center that extends for more than 4 sampling in- 
tervals (~ 2.2 arcsec) from a point near the center of the galaxy at 
twice the FWHM of the PSF, the galaxy is classified as having a 
blue core. The blue core region is defined to end where the colour 
gradients become negative. Otherwise, the galaxy is classified as 
having a red core. We fit the colour profile in the radius range from 
twice the FWHM of the PSF to where the excess blue emission 
ends (for objects with blue cores) or where the total errors reach 
0.5 mag (for objects with red cores), and define the colour gradient 
G(u - r) as the slope of the fit. Figure [6] shows the u-r profiles 
and the range of our adopted linear fits for two BCGs from the R08 
X-ray sample. 



4 RESULTS 

4.1 BCGs and control galaxies 

In FigureQ we compare the normalized distribution of colours and 
colour differences for SI BCGs (red) (see definition in Sect 12. 2t 
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Figure 6. u - r differential color profiles of two BCGs from the R08 X-ray 
sample (A1795 and A2065). The dashed lines show the 1 tr range of the 
profile. The best-fit gradients are overplotted as green straight lines in the 
fitting ranges as defined in Sect. 13. 51 The values of the best-fit gradients are 
shown in the left-down comer of each profile. 



and control galaxies (black). Colours are measured in two zones: 
inside R 50 (denoted by colour,,,) and between R 50 and R 90 (de- 
noted by colour ■„,„). Colour differences here follow the definition 
in Sect 2. The results of a Kolmogorov-Smirnov test (specifically, 
the logarithm of the K-S probability that the two distributions are 
drawn from an identical parent population) are given in each panel. 
Smaller values mean larger differences between the two distribu- 
tions (a 99 percent probability that the distributions are different 
will have a value of -2). Dashed lines in each panel indicate the 
median values of the distributions. 

The top two panels show SDSS colours measured for the 
BCGs and control galaxies. There is no systematic difference in 
median colour between the two samples. However, non-BCGs 
clearly have larger scatter in their colours than BCGs, particularly 
in their outer regions. This translates into a broader range in colour 
differences, which is consistent with the conclusions of vdL07, who 
found that BCGs have more homogeneous stellar populations than 
non-BCGs. 

There are clear systematic shifts in the NUV - r colours of 
BCGs and non-BCGs in the sense that BCGs have bluer outer 
NUV - r colours and redder inner NUV - r colours. The first 
question one might ask is whether these shifts are real, or whether 
they simply reflect that fact that BCGs have larger sizes (i?5o) and 
lower optical surface brightness than non-BCGs. As discussed in 
Sect. 13.41 this may indeed partly explain why BCGs have bluer 
outer NUV - r colours. Using the colour distributions of our sam- 
ple of BCGs and non-BCGs, we estimate that the net shift due to 
structural differences between the two kinds of galaxies is about 
0.1 mag, and thus cannot explain the differences between the two 
distributions (the reddest non-BCGs are ~0.4 mag redder than the 
BCGs). Also, we find that the inner NUV - r colours of BCGs are 
redder than those of non-BCGs, which goes in the opposite direc- 
tion to the effect predicted from the simulations In summary, our 
result that BCGs have bluer outer NUV- r colours than non-BCGs 
may be weaker than suggested in Figure|7] but the result that BCGs 
have redder inner NUV - r colours than non-BCGs is robust. This 
would imply that the inner stellar populations of BCGs are older, 
more metal rich, or more dusty (or some combination of these). 

We find no detectable differences in the FUV - NUV colours 
of BCGs and non-BCGs. This may be due to lower quality of the 
FUV images. 

We note that the BCGs and comparison galaxies have been 
matched in redshift, so ^-corrections are not an issue in our anal- 
ysis. Finally, as discussed in Sect |2.3l the control galaxies correp- 
sonding to BCGs in less massive clusters may be "contaminated" 
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Figure 7. Distribution of 2-zone colours and colour differences for optical BCGs (red) and the control sample of massive field galaxies (black). The dashed 
lines mark the median values for the two distributions. 



with BCGs. We have checked whether there is any trend in the dif- 
ferences betweeen the colours of BCGs and non-BCGs as a func- 
tion of cluster mass, but we did not find any effect. 

4.2 Optical BCGs versus X-ray BCGs 

Figure [8] shows the distributions of cr c [ uster and X- ray luminosity 
(when it is available from ROSAT. Ishen et al.l l l2008h ) for the vdL07 
sample. In the following analysis, we limit the optical BCG sam- 
ple to those objects in clusters with log (roaster > 2.8 km s" 1 . Af- 



ter including galaxies with GALEX AIS photometry, the sample 
consists of 60 BCGs. The X-ray selected BCGs from R08 (Ta- 
ble [TJ are still significantly biased towards higher cr c i us , er and even 
more severely towards higher X-ray luminosities compared to the 
optically-selected clusters. We will come back to this point later in 
our discussion. 

Figure [9] plots the BCGs at z<0.1 in a series of 
colour/magnitude and colour/colour diagrams. BCGs in optically- 
selected clusters are shown as black points and BCGs in X-ray se- 
lected clusters are shown as coloured symbols (cool-core BCGs in 
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green and non-cool-core BCGs in red). Most X-ray selected BCGs 
have high r-band luminosities compared to the optical BCGs. X- 
ray selected BCGs have similar g — r and u - r colours to the op- 
tically selected BCGs, but blue inner NUV-r colours. Among the 
X-ray selected BCGs, the cool-core BCGs are always bluer than 
the non-cool-core ones (cool-core and non-cool-core BCGs sepa- 
rate at a NUV - r colour of ~5.8). This demonstrates that the UV 
and optical colours of the X-ray selected BCGs are indeed corre- 
lated with the cooling time of the gas. One of the objects with the 
bluest NUV - r colour is t he BCG in the well studied very stron g 
cooling flow cluster A1795 dO'Dea et al.ll2004lMittaz et alj200lh . 
The other blue outlier is M87, which is an AGN. Th e synchrotron 
emiss ion from its jet is quite prominent in the UV jMadrid et al.l 
l2007h . 

In Table 2, we list K-S test probabilities for colour differences 
between optically and X-ray selected BCGs, as well as optically 
selected BCGs and cool-core BCGs. The results confirm our asser- 
tion that the most significant differences are found for UV/optical 
colours and for optically-selected and cool-core BCGs. Could the 
large differences in UV/optical colours be an artifact of our con- 
volution algorithm? As discussed in Sect 13.41 the inner NUV - r 
colours are more severely underestimated for red objects than for 
blue objects. So if we were to correct our colour measurements to 
the value before convolution, red galaxies will become redder in 
NUV - r while blue galaxies will not change much. In other words, 
the NUV - r colour differences between the cool-core and non- 
cool-core BCGs would be enhanced even more. 



4.3 Consistency with previous work 

In the previous subsection, we showed that cool-core and non-cool- 
core clusters segregated rather cleanly in colour space, implying 
that the presence of young stars in the BCG is related to the cooling 
time of the gas in cluster cores. 

There have been a number of studies of X-ray selected clus- 
ters that have found that there is a central entropy (or l coo i) thresh- 
old be low which BCGs have significant Ha line emission and blue 
cores teildfell et alj|2007l ; iRaffertv et ai]|2008l ; IPipino et alj|2008l ; 
ICavagnolo et alj2008al) . 

We now check whether we reproduce these results using the 
sample of X-ray BCGs from R08. For this analysis, we do not im- 
pose a cut on the sample at z <0. 1 - the reason is that this would 
reduce the sample size by too much. Table [Tj shows the values of 
the u — r colour gradients given in R08 and our measurements 
from the SDSS images. The measurements are generally consis- 
tent given the uncertainties . When available, we t ake the Ha fluxes 
from Crawford et al.l J 19991) (to be consistent with lCavagnolo et"al] 
J2008al) ). otherwise, we adopt the Ho- fluxes in the SDSS DR7 
emission line catalog released by the MPA/JHU group. 

We plot it - r colour and Ho- luminosity versus t„, / in Fig- 
ures[lO]and 1 1 1 [respectively. We use the cross correlation parameter 
p = {(X - X)(Y - ?)) to evaluate the significance of corrleation be- 
tween parameter X and Y. As can be seen, we do verify a threshold 
t C ooi value of ~ 1 Gyr below which some BCGs exhibit significant 
Ha luminosity and blue cores. It should be noted, however, that not 
all BCGs with t co „i < 1 Gyr are strongly star-forming. There is a 
significant scatter in Ha luminosity and colour gradient values and 
only weak correlation with the actual value of t cool below the thresh- 
hold value. If we plot u — r colour gradient versus Ho- luminosity 
for these X-ray BCGs (Figure[T2l. there is a significant correlation 
between the two quantities (albeit with large scatter). We also find 
similar correlations between u — r colour differences (following the 
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Figure 10. The u — r colour gradient as denned in R08 (bottom, see text) is 
plotted as a function of the cooling time of the gas in units of 10 s years. The 
dashed line marks 0.8 Gyr, which R08 claim to be the threshold for finding 
blue cores. Numbers are used to mark individual BCGs as listed in Table 1 . 
p is the cross correlation parameter, as denned in Section 4.3. . 



definition from R08) and both Hey luminosity and Ho- equivalent 
width. 

We note that the scatter may arise because Ho- line emission is 
produced by ionizing photons from massive stars with lifetimes of 
~ 10 7 years, whereas the u - r colour gradients will be sensitive to 
star formation that has taken place over Gyr timescales. In addition, 
it is possible that part of the Ho- line emission does not originate 
from HII regions, but is excited by a central AGN. 

We have also checked whether there is a correlation between 
u — r colour gradient and Ha luminosity in the ROSAT (top) and 
GALEX MIS (bottom) detected sample of optically-selected BCGs 
from vdL07 and the results are shown in Figure [T3l the lower red- 
shift range of z < 0.1 ensures that the u - r colour gradients can 
be measured relatively accurately. In the top panel, the points are 
colour-coded according to the total X-ray luminosity of the clus- 
ter (clusters with ROSAT X-ray luminosities greater than 10 43 5 erg 
s -1 are plotted in red). In the bottom panel, the points are coded 
according to the NUV - r colour of the BCG. We fail to find any 
correlation between colour gradients and Ha luminosity for any of 
the subsamples defined from the optically-selected clusters. 

This may not be too surprising in view of the fact that the op- 
tically selected BCGs span a much smaller range in Ha luminosity, 
(most L Ha <l.lxl0 40 erg/s). Much of the correlation seen in Fig- 
ure[l2]is driven by the BCGs with the very higher Ha luminosities. 
Over the range in Ha lumnosity where the R08 sample overlaps 
the optically selected clusters, the results are actually in fairly good 
agreement. 



5 SUMMARY AND DISCUSSION 

The main aim of this paper is to investigate the claim that the bright 
galaxy located at the center of a massive dark matter halo (BCG) 
forms stars at a higher rate than galaxies of similar mass and struc- 
ture that are not located at cluster or group centres. Studies of BCGs 
from X-ray selected cluster samples have generally supported this 
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Table 1. Properties of the BCGs in the X-ray selected sample of (R08). Columns 2 and 3 list the velocity dispersions and X-ray luminosities of the host 
cluster s, which ha ve been obtained fro m a variety of different sources as indicated by a superscript to the values of L x as follows: s r efers to lsmith et al] 
1200 ll). e refers tolEbeling et all ll996t). es re f ers tolEgami et alj 120061). h refers tolHashimoto et all d2007n . H refers to lHoekstra et alj <2007l) . HI refers to 
iHelsdon et al] J200lh . p refers to lPeres e t alj 1 19 981). w refer s to lD. A.White et al" Hl997l) . Column 4 lists the Ho- emission line luminosities of the BCGs. c 
means the Ho- luminosity comes from lCrawford et alj Jl999t) . otherwise it is from SDSS DR4. Column 5 is the cooling time of the gas calculated by R08. 
Column 6 lists the u - r colour gradient given by R08, while Column 7 is our own measurement of the same quantity using imaging data from SDSS 




Figure 8. cr c i uster (right panel) and X-ray luminosity (left panel) of the entire optical BCG sample (black) and subset of the optically-selected BCGs with data 
from ROSAT (blue). The diamonds mark the values of the X-ray selected R08 BCGs used in this work. 
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Figure 9. Comparison between the inner and outer colours of optical BCGs (black, filled circles; only those in clusters with o"„ >2.8 are plotted), cool-core 
(green, filled triangle) and non-cool-core BCGs (red, crosses). 
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Table 2. K-S test values (the logarithm of the K-S probability that the two distributions are drawn from an identical parent population) for differences between 
optically selected BCGs, X-ray selected BCGs and cool-core X-ray BCGs. It shows that the difference of (NUV - r),„ between X-ray selected (especially 
cool-core) and optical selected BCGs is strong. 



claim. Anlayses of BCGs from optically-selected cluster samples 
have failed to find any significant enhancement in star formation 
when compared to control samples of non-BCGs. 

Recent work on optically-selected BCGs relied on spectro- 
scopic indicators derived from Sloan Digital Sky Survey spectra. 
This has the disadvantage of only probing the most recent star for- 
mation that is occurring in the cores of the BCGs. In this paper, we 
analyzed a sample of 1 13 BCGs with z < 0. 1 with optical imaging 
from SDSS and UV data from GALEX. We convolved the SDSS 
images to match the resolution of the GALEX data and measured 
UV/optical colours both in the inner and outer regions of the BCGs. 
We did the same for a smaller sample of X-ray selected BCGs with 
gas cooling times derived from Chandra data. 

Our main conclusions concerning our sample of optically- 
selected BCGs are the following: 

(i) We confirm the conclusion of vdL07 that optically-selected 
BCGs have similar median g-r and u—r colours as control samples 
of non-BCGs with the same stellar mass and concentration. This 
conclusion holds for colours measured in both the inner and the 
outer regions of the galaxies. 

(ii) Optically-selected BCGs have smaller scatter in their g-r 
and u — r colours than the control non-BCGs. 



(iii) Optically-selected BCGs have slighly redder NUV - r 
colours in their centres than the control galaxies. 

(iv) Correlations between G(u - r) and Ha luminosity are very 
weak. 

Taken together, the three findings lead us to conclude that optically- 
selected BCGs are slightly older and have had less recent star for- 
mation than the galaxies in the control sample. 

Our main conclusions concerning our sample of X-ray BCGs 
are the following: 

(i) The optical (g - r and u - r) colours of the BCGs in the X-ray 
selected cluster sample span a similar range as those in the subset of 
optical clusters with high cluster (log cr > 2.8) velocity dispersions. 
Their NUV - r colours, on the other hand, are significantly bluer. 

(ii) Among the X-ray selected BCGs, those that are located in 
the centre of a cool-core cluster (defined to have central cooling 
times of less than 1 Gyr) always have bluer optical and UV/optical 
colours than those that are located in clusters where the central gas 
cooling times are long. 

(iii) We confirm that there appears to be a threshold value of the 
central cooling time below which star formation does occur in a 
subset of clusters. 

Our main conclusion, therefore, is that the location of a galaxy 
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Figure 11. Ha line luminosity is plotted against t COB \ for the R08 BCG sam- 
ple. Numbers are used to mark individual BCGs. p is the cross correlation 
parameter. 
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Figure 12. The correlation between G(u - r) and i/cr line luminosities for 
the R08 sample. Numbers are used to mark individual BCGs. p is the cross 
correlation parameter. 



at the center of a dark matter halo is not sufficient to determine 
whether or not it is currently forming stars - one must also have 
information about the thermo-dynamic state of the gas in the core of 
the dark mat t er halo . These results agree with other similar studies. 
lEgami et al.1 Eo0& found that the majority of the BCGs are not 
particularly in frared-luminous com pared with other massive early- 
type galaxies. lEdward et all J2007h found their red BCGs selected 
from the SDSS sample had a low fraction with emission lines (1 1±2 
percent), comparable to that found in their control sample However, 
both papers noted that for BCGs in clusters where cooling times 
are very short and the predicted cooling flow rates are high, star 
formation is prominently enhanced. 
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Figure 13. The correlation between G(« - r) and Ha line luminosities for 
the optically selected BCG of vdL07. In the bottom panel, the BCGs with 
NUV - r <5.8 are in green (filled triangles), and others are in red (crosses). 
In the top panel, the BCGs with ROSAT X-ray luminosities > 10 43 ' 5 erg 
s are in red (pluses), otherwise they are in grey (solid points). Blue points 
(diamonds) indicate R08 BCGs. p is the cross correlation parameter for the 
symbols of corresponding colour. 



Perhaps the most important lesson to take away from this anal- 
ysis, is that BCGs selected from optical cluster surveys and BCGs 
selected from X-ray cluster surveys will be different. Low redshift 
spectroscopic surveys such as the SDSS are designed to provide a 
complete census of the most massive galaxies out to redshifts of 
0.2 - 0.3. The fact that only a relatively small percentage of BCGs 
in these surveys have signs of recent star formation must indicate 
that the majority of nearby BCGs do not sit at the centers of dark 
matter halos in wh ich central coolin g time of the gas is less than 1 
Gyr. As shown bv lBest et al] J2007I) . a large fraction (20-30%) of 
optically-selected BCGs host radio-loud AGN, which may play a 
key role in preventing the gas from cooling. 

It remains to be seen whether the same conclusion will hold up 
at higher redshifts. Many of the most strongly star-forming BCGs 
in the current X-ray samples are located at relatively high redshifts. 
We checked the full R08 sample and found that the ratio of blue 
core BCGs to those with negative gradients increases with redshift: 
19:27 for the whole sample, 13:8 for BCGs at z >0 .1, and 6:3 for 
BCGs at z >0.2. In addition. ICrawford et al.1 jl999h found that the 
fraction of emission line BCGs increases with redshift. We note 
that the cluster sample of vdL07 is defined to lie at redshifts below 
0.1, so this may explain why the BCGs in this sample are predom- 
inantly inactive. On the other hand, only the most X-ray luminous 
clusters can currently be detected at higher redshifts and these are 
likely to be biased to systems with short central cooling times. A 
new generation of both optical spectroscopic surveys (e.g. SDSS- 
III's Baryon Oscillation Spectroscopic Survey (BOSS )) and deeper 
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X-ray imaging surveys over large areas of the sky (e.g. the extended 
Roentgen Survey with an Imaging Telescope Array (eROSITA) 
will be needed before these issues can be fully understood. 
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